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i—i • V407 Cyg was, before 2010 Mar., known only as a D-type symbiotic binary system in which the Mira variable has a pulsation 

j/j ' period of approximately 750 days, one of the longest known. On 2010 Mar 10, it was discovered in outburst, eventually reaching 

Ci , V < 8. This is the first recorded nova event for this system, but it closely resembles the spectroscopic development of RS Oph, the 

prototypical symbiotic-like recurrent nova. It was also detected by Fermi above 100 MeV and displayed strong, likely nonthermal 
centimeter wavelength radio emission. Unlike classical novae occurring in compact cataclysmic binary systems, for which the ejecta 
undergo free ballistic expansion, this explosion occurred within the dense, complex wind of a Mira variable companion. This paper 



^ , concentrates on the development of the shock and its passage through the Mira wind. We also present some constraints on the binary 

system properties. Using medium and high resolution ground-based optical spectra, visual and Swift UV photometry, and Swift 
\l ' X-ray spectrophotometry, we describe the behavior of the high-velocity profile evolution for this nova during its first three months. 

' Using the diffuse interstellar bands visible in the high-resolution optical spectra, we obtain an extinction E(B - V) s;0.45±0.05. The 

, spectral type of the red giant during this period, when the star was at R minimum, was no earlier than M7 III. The peak of the X-ray 

. emission occurred at about day 40 with a broad maximum and decline after day 50. The main changes in the optical spectrum began 

at around that time. The He II 4686A line first appeared between days 7 and 14 and initially displayed a broad, symmetric profile that 
is characteristic of all species before day 60. The profile development thereafter depended on ionization state. Low-excitation lines 
remained comparatively narrow, with Viini max of order 200-400 km s"'. They were systematically more symmetric than lines such as 
[Ca V], [Fe VII], [Fe X], and He II, all of which showed a sequence of profile changes going from symmetric to a blue wing similar 
to that of the low ionization species but with a red wing extended to as high as 600 km s"'. The [O I] 6300, 6364 doublet showed 
a narrow wind-emission component near the rest velocity of the system and a broad component, 200-300 km s"', whose relative 
, ; , intensity increased in time. Forbidden lines of N II and O III had two separate contributors to the profiles, a broad line that increased 

■ in strength and velocity width, exceeding 200 km s"', and narrow components from a surrounding ionized region at higher velocity 

than the Mira wind. The Na I D doublet developed a broad component with similar velocity width to the other low-ionization species. 
The O VI Raman features observed in recent outbursts of RS Oph were not detected. We interpret these variations as aspherical 
expansion of the ejecta within the Mira wind. The blue side is from the shock penetrating into the wind while the red wing is from the 
low-density periphery. The maximum radial velocities obey power laws, Vn,ax ~ f " with na:l/3 for red wing and ~0.8 for the blue. 

Key words. Stars-individual(V407 Cyg, RS Oph), symbiotic stars, physical processes, novae 

1 . Introduction Darnley 2008; Bode 2010). Both show repeated episodes of ex- 

plosive mass ejection, between lO^^M© and lO^^M©, on recur- 
The recurrent novae (RNe), the rarest and most extreme of the '"e"'^^ timescales of decades with velocities (based on the ultra- 
classical novae, ai-e known in two distinct varieties (for a gen- "'^^^^ resonance transitions) that can exceed 5000 km s"'. One 
eral series of overviews and specific details on observations and g™up is indistinguishable from classical novae: compact sys- 
models see Bode & Evans 2008; Evans, Bode, O'Brien & M.J. ^ems with orbital period <1 day and a low-mass (normally mam- 
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sequence) star that transfers mass to its companion white dwarf 
(WD) by Roche lobe overflow in a stream and accretion disk 
(Shaefer 2010). The other even rarer group, far less well stud- 
ied but no less significant, resembles the symbiotic binaries, the 
so-called "symbiotic-like recurrent novae" (SyRNe), in which 
the mass-losing star is a red giant (RG) (e.g. Gonzalez-Riestra 
1992; Anupama & Sethi 1994; Shore et al. 1996; Shore 2008; 
Anupama 2008; Schaefer 2010). Until this year, there were only 
four known in the Galaxy: V745 Sco, V3890 Sgr, RS Oph, and 
T CrB. Unlike the compact systems (e.g. U Sco), which have 
orbital periods coincident with classical cataclysmics, these sys- 
tems have long orbital periods (years) and the WD likely ac- 
cretes mass from the RG wind, instead of a stream. T CrB and 
RS Oph, for which accurate orbital parameters are known, host 
a WD whose mass exceeds that of the RG and for which the de- 
generate is more massive than in the symbiotic novae (not to be 
confused, these do not explosively eject mass and, instead, dis- 
play a thermonuclear runaway that enters a long-term period of 
burning with an inflated WD envelope). The latter type are on 
the more active end of the symbiotic systems, the WD is the less 
massive component and <1 M©. 

Recurrent novae are produced by the initiation of a ther- 
monuclear runaway provoked by accretion on massive WDs, 
which are nearer to the Chandrasekhar limit (for RS Oph, 
MwD ~ 1.2 M0) than for normal symbiotics or symbiotic novae 
(< 1 M0) (see Stan-field 2008; Anupama 2008, Mikolajewska 
2008, Brandi et al. 2009). The connection of these systems, and 
their supersoft XR emission after outburst, to SN la progenitors 
has long been a scenario (e.g. Munari & Renzini 1992; Hernanz 
& Jose 2008; Sokoloski et al. 2006; Walder et al. 2008; Bode 
2010; Schaefer 2010). The amount of mass that remains on the 
WD after the explosion is the critical determining factor. It is still 
unknown if there is a wind subsequent to the explosion that con- 
tinues to remove the accreted mass, thus preventing the eventual 
critical accumulation to produce a collapse of the WD, or if the 
explosion ends abruptly and ejects only a fraction of the mass. 
However, no symbiotic had ever been observed to display a clas- 
sical nova explosion until that of V407 Cyg. This is why the 
system is both unique and significant. 

V407 Cyg, a D-type symbiotic binary with a long period 
Mira RG (Munari, Margoni, & Stagni 1990), underwent its first 
ever recorded explosive mass ejection on 2010 Mar 10, reaching 
an unfiltered magnitude of 6.8 (Nishiyama & Kabashima 2010) 
and V magnitude peak of 7.9. Historically, V407 Cyg was first 
identified by its nova-like behavior in 1936 (Hoffmeister 1949, 
Meinunger 1966) in the nipg range 14 to 16.5, but with identified 
contributions to the light curve from the Mira variable. Since 
that time, it has been observed to periodically increase in bright- 
ness according to the long-period Mira, and minor, slowly de- 
veloping outbursts, never becoming brighter than approximately 
y = 1 1.5 during the years 1984 - 1995 (Kolotilov et al. 1998). A 
slow-developing outburst in the late 1990s was identified from 
the light curve (Kolotilov et al. 2003, Shugarov et al. 2007), as 
it was claimed to become brighter than the outburst of 1936, but 
never being brighter than V ~ 11.8. Large gaps in the light curve 
exist from the discovery data epoch until around 1990, which are 
of sufficient duration to have missed any outburst similar to the 
current event. Unfiltered optical monitoring of V407 Cyg dur- 
ing the two years prior to the March outburst shows a step-like 
brightness increase of nearly three magnitudes over a one year 
interval prior to outburst (Abdo et al. 2010). The 2010 outburst 
increased the brightness by an additional 2.5 magnitudes. These 
step-like brightness increases, along with other irregular bright- 
ness variations, are vaguely present in the historical light curve. 



Spectra from 2010 March 14-16, reported by C. Buifl 
showed very broad emission lines for all permitted transitions, 
especially the H Balmer lines, and He I, with an ejection velocity 
>2000 km s"', with superimposed very narrow emission lines. 
The surprise that V407 Cyg was previously considered a classi- 
cal - although somewhat unusual - symbiotic was compounded 
by the high-energy detection of both X-rays (XR) and y-rays 
above 100 MeV (Abdo et al. 2010). These were interpreted as 
the signature of 7r"-production by shock-energized protons in the 
RG wind during the first four days of the outburst, when the ve- 
locities exceeded 2000 km s"' . 

In this first paper we focus on the shock and its development. 
In subsequent papers we will present the analysis of the circum- 
stellar material illuminated by the photoionization produced dur- 
ing the outburst and a line list constructed using high resolution 
spectra obtained from the Nordic Optical Telescope (NOT) on 
day 22 of the event, before the peak of the XR emission. 



2. Observations 

Our optical observational data set consists of spectra taken 
between 2010 Mar 24 and 2010 Jul 16 with the Ondfejov 
Observatory Zeiss 2.0 m telescope coude spectrograph and with 
the 2.6 m Nordic Optical Telescope (NOT) fiber-optic echelle 
spectrograph (FIES, program P40-423). The Ondfejov spectra 
taken with the SITe005 800x2000 chip were mainly obtained 
in the vicinity of Ha, supplemented by spectra at both bluer 
and redder wavelengths. The dispersion of 0.24A per pixel pro- 
vided spectral coverage of approximately 500 A, and exposure 
times ranged from 60 s to 6700 s. The FIES spectra were ob- 
tained with a dispersion of 0.023A px"' in high-resolution mode 
covering the spectral interval from 3635 to 7364A and 0.035A 
px"' in medium-resolution mode covering the spectral interval 
from 3680 to 7300A. Exposures ranged from 100 s to 3000 s. 
Absolute fluxes for the NOT spectra were obtained using the 
flux standard star BD -h28°421 1 observed on Jun 3 and Jun 23 at 
high resolution. Other spectra were not absolutely calibrated. All 
NOT spectra were reduced using IRAF, FIESTools, and our own 
special purpose routines written in IDlQ In several instances, 
contemporaneous spectra between the two observatories allow 
us to correlate and cross-calibrate the data. The journal of our 
observations is presented in Tables la and lb, which provides 
the observation start time in the systems Universal Time, Julian 
Date, JD, and AT, the approximate time since maximum optical 
light (TO). The TO parameter has been adopted to be 2010 Mar 
10.813 = JD 2455266.313 (see section 3.1) forXR analysis. For 
the journal of observations we adopt the value JD 2455266.30 
as a guide to following the discussion, based on the three nearly 
coincident reportings of maximum brightness in CBET #2199. 

Note: (a) TO (JD) = 2455266.30; (b) HR = high resolution, 
MR = moderate resolution Note: (a) AT = T - TO, with TO (JD) 
= 2455266.30 



' URL: http://www.astrosurf.coin/~buil/v407cyg/obs.html 
- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
US National Science Foundation. 
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Table 1. Table la. Nordic Optical Telescope Observations 



Spcclrum ID 


Dale 




UT 






JD-240000 


AT (days)" 


Uxp (sec) 


Object" 


FItc310098 


2010-04- 


-01 


05: 


:28: 


:24. 


,8 


55287.72806 


21.4 


1000 


V407 Cyg HR 


FItd010076 


2010-04- 


-02 


05: 


:38: 


:54. 


.7 


55288.73536 


22.4 


125 


V407 Cyg HR 


FItdO 10077 


2010-04- 


■02 


05: 


:43: 


:11, 


,5 


55288.73833 


22.4 


400 


V407 Cyg HR 


FItd230110 


2010-04- 


■24 


05: 


:11 


:12, 


,5 


55310.71612 


44.4 


1000 


V407 Cyg HR 


FItd270157 


2010-04- 


■28 


05: 


:30: 


:37, 


,4 


55314.72960 


48.4 


100 


V407 Cyg HR 


FItd270158 


2010-04- 


■28 


05: 


:33: 


:40, 


,3 


55314.73172 


48.4 


1000 


V407 Cyg HR 


FIte 170040 


2010-05- 


■18 


04: 


:52: 


:22, 


,9 


55334.70304 


68.4 


4 


V407 Cyg HR 


FItel70041 


2010-05- 


■18 


04: 


:53: 


:48. 


,6 


55334.70403 


68.4 


100 


V407 Cyg HR 


FItel70042 


2010-05- 


-18 


04: 


:56: 


:54. 


,0 


55334.70615 


68.4 


1000 


V407 Cyg HR 


FIte260044 


2010-05- 


-27 


04: 


;38: 


:28. 


,8 


55343.69339 


77.4 


400 


V407 Cyg HR 


FIte260045 


2010-05- 


-27 


04: 


:46: 


:30. 


.5 


55343.69896 


77.4 


3000 


V407 Cyg HR 


FItf020080 


2010-06- 


■03 


04: 


:03: 


:20, 


,1 


55350.66898 


84.4 


200 


V407 Cyg HR 


FItf020081 


2010-06- 


■03 


04: 


:08: 


:01, 


,8 


55350.67224 


84.4 


2564 


V407 Cyg HR 


FItf020084 


2010-06- 


■03 


04: 


:57: 


:26, 


,2 


55350.70655 


84.4 


500 


BD+28d4211 HR 


FItf 150084 


2010-06- 


■16 


04: 


:08: 


■29. 


,0 


55363.67256 


97.4 


200 


V407 Cyg HR 


Fltf 150085 


2010-06- 


■16 


04: 


:13: 


■Al. 


,8 


55363.67583 


97.4 


3000 


V407 Cyg HR 


FItf220081 


2010-06- 


■23 


03: 


:10: 


■23. 


,3 


55370.63221 


104.3 


500 


BD-l-28d4211 MR 


FItf220083 


2010-06- 


-23 


03: 


:27: 


:37. 


,6 


55370.64419 


104.3 


270 


V407 Cyg MR 


FItf220084 


2010-06- 


-23 


03: 


:33: 


:29. 


,5 


55370.64826 


104.3 


2700 


V407 Cyg MR 


FItgl50076 


2010-07- 


-16 


00: 


;00: 


:52. 


,9 


55393.50061 


127.2 


270 


V407 Cyg MR 


FItgl50077 


2010-07- 


-16 


00: 


:06: 


:44, 


,7 


55393.50468 


127.2 


2700 


V407 Cyg MR 



Table 2. Table lb. Ondfejov Observations 



Spectrum ID 


Date 




UT 


JD - 2400000. 


AT (days)" 


texp (sec) 


^min ( 


tc230055 


2010- 


■03^ 


-24 


01: 


: 45:46 


55279, 


,57345 


13, 


.3 


1641 


6258 


tc230056 


2010- 


■03^ 


-24 


02: 


:13:45 


55279, 


,59288 


13, 


.3 


373 


6258 


tc230058 


2010- 


■03^ 


-24 


02: 


:24:21 


55279, 


,60024 


13, 


.3 


447 


6258 


tc230059 


2010- 


-03- 


-24 


02: 


:32:24 


55279. 


,60583 


13, 


.3 


389 


6258 


tc230060 


2010- 


-03- 


-24 


02: 


:39:31 


55279. 


,61078 


13, 


.3 


360 


6258 


tc230065 


2010- 


-03- 


-24 


03: 


:07:52 


55279. 


,63046 


13, 


.3 


328 


8392 


tc230066 


2010- 


-03- 


-24 


03: 


: 14:45 


55279. 


,63524 


13 


.3 


1985 


8392 


td020076 


20 10^ 


■04^ 


-03 


01: 


:59:53 


55289, 


,58325 


23, 


.3 


577 


6258 


td020077 


20 10^ 


■04^ 


-03 


02: 


: 10:07 


55289, 


,59036 


23, 


.3 


586 


6258 


td020078 


20 10^ 


■04^ 


-03 


02: 


: 20:32 


55289, 


,59759 


23, 


.3 


599 


6258 


td030034 


20 10^ 


■04^ 


-04 


01: 


:33:58 


55290, 


,56525 


24, 


.3 


541 


6257 


td030035 


20 10^ 


■04^ 


-04 


01: 


:43:38 


55290, 


,57197 


24, 


.3 


586 


6257 


td030036 


2010- 


-04- 


-04 


01: 


:54:04 


55290. 


,57921 


24, 


.3 


626 


6257 


td050019 


2010- 


-04- 


-06 


00: 


:35:32 


55292. 


,52468 


26, 


.2 


936 


6257 


td050022 


2010- 


-04- 


-06 


00: 


:58:26 


55292. 


,54058 


26, 


.2 


3759 


4753 


td060048 


2010- 


-04- 


-06 


23: 


: 16:59 


55293. 


,47013 


27, 


.2 


1255 


6257 


td060055 


2010- 


■04^ 


-06 


23: 


:48:47 


55293, 


,49221 


27, 


.2 


6718 


4797 


td060063 


2010- 


■04^ 


-07 


01: 


:51:17 


55293, 


,57728 


27, 


.3 


1996 


6400 


td070033 


2010- 


■04^ 


-08 


02: 


:33:25 


55294, 


,60654 


28, 


.3 


60 


6257 


td070034 


2010- 


■04^ 


-08 


02: 


:35:58 


55294, 


,60831 


28, 


.3 


300 


6257 


td070037 


2010- 


■04^ 


-08 


02: 


:48:24 


55294, 


,61694 


28, 


.3 


60 


8392 


td070038 


2010- 


-04- 


-08 


02: 


:50:52 


55294. 


,61866 


28, 


.3 


2000 


8392 


td210056 


2010- 


-04- 


-22 


02: 


:37:40 


55308. 


,60949 


42, 


.3 


500 


8395 


td210057 


2010- 


-04- 


-22 


02: 


:47:17 


55308. 


,61617 


42, 


.3 


500 


6259 


td260020 


2010- 


-04- 


-26 


23: 


: 12:42 


55313. 


,46715 


47, 


.2 


4135 


6261 


td260021 


2010- 


■04^ 


-27 


00: 


:25:59 


55313, 


,51804 


47, 


.2 


1106 


6261 


td280030 


2010- 


■04^ 


-28 


23: 


:37:32 


55315, 


,48440 


49, 


.2 


800 


6260 


td280036 


2010- 


■04^ 


-29 


02: 


:02:32 


55315, 


,58509 


49, 


.3 


1999 


8396 


tf040013 


2010- 


■06^ 


-04 


23: 


:06:20 


55352, 


,46273 


86, 


.2 


3000 


6258 


tf040020 


20 10^ 


■06^ 


-05 


00: 


:15:17 


55352, 


,51061 


86, 


.2 


300 


8394 


tf040022 


2010- 


-06- 


-05 


00: 


:25:57 


55352. 


,51802 


86, 


.2 


2400 


8394 


tf040030 


2010- 


-06- 


-05 


01: 


: 19:26 


55352. 


,55516 


86, 


.3 


2000 


7681 


tf060016 


2010- 


-06- 


-06 


20: 


:46:29 


55354, 


,36561 


88, 


.1 


4499 


6260 


tf070017 


2010- 


-06- 


-07 


20: 


:35:33 


55355. 


,35802 


89 


.1 


3000 


6740 


tf070021 


20 10^ 


■06^ 


-07 


21: 


:32:43 


55355, 


,39772 


89, 


.1 


3600 


7246 


tf070030 


20 10^ 


■06^ 


-07 


23: 


:31:08 


55355, 


,47995 


89, 


.2 


2100 


6258 


tfl40016 


2010- 


■06- 


-14 


21: 


:24:24 


55362, 


,39194 


96 


.1 


3600 


8393 


tfl40020 


2010^ 


■06^ 


-14 


22: 


:32:02 


55362. 


,43891 


96 


.1 


1800 


6256 


tgl20015 


2010- 


-07- 


-12 


21: 


: 16:22 


55390, 


,38637 


124.1 


3600 


6259 
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3. Analysis 

3.1. Swift Observations: XRTand UVOT 

Swift (see Gehrels et al., 2004) started observing V407 Cyg three 
days after the peak of the optical light-curve (2010 Mar 10.813 
= JD 2455266.313, hereafter denoted as TO), with both the X- 
ray Telescope (XRT, Burrows et al., 2005) and the UV/Optical 
Telescope (UVOT, Roming et al., 2005). Short observations (1 
to 2 ksec) were performed every few days until Apr 2 (TO-i-23.1 
days, JD 2455280, the start of our Ondfejov observations), at 
which point the exposure time was increased to ~4 ksec every 
other day (with a ~10 ks observation on 2010 Apr 17; TO-l-37.6 
days). From 2010 Apr 28 until May 6 (days 48.6-56.3), 3 ksec 
were obtained every four days. From 2010 May 18 (TO-i-68.7 
days), observations of about 1 ksec were taken weekly until 2010 
July 2.7. Data were processed using version 3.5 of the Swift soft- 
ware (released as part of HEASoft 6.8 on 2009 December 03) 
and the most up-to-date calibration files. 

Initially the XRT PC mode count rate remained low, ~0.01- 
0.02 count s"^ until TO-i-15.2 days. The rate started to increase at 
T0-H12.4 days, reaching a peak of ~0.26 count s'^ on TO-i-30.5 
days. The X-rays then began a slow decay until around TO-i-56 
days, at which point the count rate stopped decreasing, remain- 
ing between 0.08-0.09 count s ' . The light curve is shown in Fig 
1. 

The UV source was bright and most of the data were col- 
lected in the uvm2 filter (central wavelength 2246A); observa- 
tions in u and uvwl were saturated and are not shown. The ini- 
tial data showed a slow decline (approximately 0.02 mag day"^) 
that steepened to ~0.09 mag day^' around 40 days after peak; af- 
ter TO+56 days the decay once more slowed to 0.05 mag day"'. 
This reduced rate of decline in the U V light curve starts at around 
the same time as the XR flux stops declining, a similar feature 
appears in the AAVSO V band light curve at this time. As we 
discuss below, this is the same behavior derived from the optical 
emission lines and, in particular, from the third light continuum 
detected using the photospheric absorption lines from the Mira. 

The XR spectrum remained relatively hard throughout the 
observations, with no strong emergence of super-soft emis- 
sion. The hardness ratio (defined as the count rate ratio [1-10 
keV]/[0.3-l keV]) appears to soften until around day 30, af- 
ter which time it remains approximately constant (see Fig 1). 
Spectra were extracted before TO-i-30 days (as the XR Ught- 
curve was brightening), T0=30-60 days (as the X-rays were fad- 
ing) and after day 60 (when the XR count rate was approxi- 
mately constant). The earliest two spectra are very flat and can- 
not be fit with simple models; by fitting we ruled out all 
solar abundance models from single-temperature optically-thin 
models with a single absorption component to two-temperature 
models with two absorption components, see Table 2. We used 
the Wilms et al. (2000) abundances and absorption model and 
the mekal plasma emission model. In contrast, non-solar abun- 
dance two-temperature models fit these two spectra well. For 
the late-time (after day 60) spectrum, orily a single tempera- 
ture component was needed; however, a non-solar abundance 
was still a significant improvement. We were able to achieve 
a good fit to all three spectra simultaneously, using the same 
abundances for each, but with other parameters taking different 
values. While two emission components were required for the 
spectra before day 60, we could not distinguish between a black- 
body and an optically-thin plasma model for the lower tempera- 
ture component. Both spectral components were absorbed by a 
single-component cold-absorption model, with the colurmi den- 




se 100 
day since optical peak (201 March 1 0.81 3) 



Fig. 1. Swift XRT 0.3-10 keV count rate, hardness ratio, and 
UVOT uvm2 (2246A) magnitudes of V407 Cyg as a function 
of time since the optical peak. For comparison with the optical 
spectroscopy, note that TO = MJD 55266. 



sity decreasing with time. The temperature of the hotter (kTsi2- 
4 keV) component was also lower in the second spectrum, al- 
though the luminosity was unchanged. The (single) temperature 
of the spectrum extracted after day 60, when the XR count rate 
remained close to constant, was again lower, while the absorbing 
column was consistent with the day 30-60 spectrum. The lumi- 
nosity of the low temperature component is not well constrained 
due to the combination of its low temperature and high absorb- 
ing column. The abundances implied by this fit are very high 
for N and O, abundances for elements other than those hsted 
were fixed at solar. The simultaneous spectral fit had ;t'^/dof = 
480/383 (reduced - 1.25) for the optically-thin plus black- 
body model, and x^/dof = 489/385 (reduced ^ j 27) for die 
two-temperature opticaUy-thin model. The results of the former 
fit are shown in Table 3 and Fig 2; aU errors and limits are at 
90% confidence. The lower temperature component is respon- 
sible for 14% and 11% of the counts in the early and middle 
spectra, respectively. It remains possible that a yet more com- 
plex temperature distribution would reduce the implied high N 
abundance, although we were not able to find a sensible better 
fit with three temperature components. Consistent with the XRT 
spectral fits there is no detectable signal in the Swift BAT, with 
the 15 to 50 keV count rate remaining below 0.003 count cm~^ 
s-i. 

For both spectra, the luminosity at 2.7 kpc (Munari et al. 
1990) is approximately 1.2 x 10^"* erg s"' for days 0-30 and 30- 
60. The day 60-i- spectrum has a corresponding luminosity of 
about half this value, 6.4 x 10^^ erg s"'. This is the unabsorbed 
luminosity, for the hotter component* only. The final non-solar 
abundance enhancements tied between both spectra, with respect 
to Wilms et al. (2000) are N > 251, O = 35+^°, Ne <3.1, Na = 
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Fig. 2. The joint best fit to the Swift XRT spectra accumulated 
before day 30, from day 30 to day 60, and after day 60 from 
TO. The model is a strongly absorbed optically thin high temper- 
ature component with strongly enhanced nitrogen and oxygen 
together with a lower temperature blackbody for the first two 
spectra. See text for details. 



103^^4^, and Mg < 2.3. The upper limit to the Ne and uncertain 
O abundances do not indicate an ONe white dwarf as the site of 
the explosion. They do not, however, exclude it since the ejecta 
rapidly mix with the environmental material (Walder et al. 2008) 
the abundances of which are closer to solar than expected from 
the nucleosynthesis occurring during the nova event. Thus, the 
final elemental mixture of the ejecta may be substantially altered 
during their passage through the Mira wind. 

Table 2: x^/'iof values for solar abundance models for the XRT 



spectra 


Model 


< day 30 


day 30-60 


> day 60 


mekal*tbabs 
mekal*pcfabs 
(mekal-i-mekal)*tbabs 
(mekal-i-mekal)*pcfabs 


514/129 
263/128 
210/127 
208/126 


622/246 
615/245 
476/244 
301/243 


580/244 
580/243 
430/242 
414/241 



Table 3: Best fit non-solar abundance models for the XRT 
spectra 



spectral type of M6 III, Munari et al. (1990) derived an extinc- 
tion of E(B - V) - 0.57. We have instead used optical diffuse 
interstellar band (DIB) identifications fromHerbig (1975, 1988) 
and Jenniskens & Desert (1994). We used the calibrations for 
the DIBs given in the latter reference. Several are either blended 
with emission lines (e.g. 4428A) or too weak to be measured, 
but most of the important features are accessible. Our measure- 
ments are listed in Table 4, the uncertainty is +5mA. The re- 
sulting reddening is E(B - y)=0.45±0.09 that is consistent with 
but lower than the value derived from the energy distribution. 
We note, however, that the Mira was at minimum light during 
this outburst. We have compared the NOT spectra after MJD 
55370 with those of giants in the range M6 - M8 taken from the 
ESO/VLT/UVES Library of High Resolution Spectra of Stars 
Across the H-R Diagram archived data set (Bagnulo et al. 2003) 
and find that the present spectrum is a far better match to a 
later spectral type, close to M7 III or M8 III. In light of the ap- 
proximate agreement of these different methods, we will adopt 
E(B - y)=0.5 with a range of ±0.05. 

Table 4: DIBs in the NOT spectrum of V407 Cyg, MJD 55288 



3.2. Interstellar extinction and spectral type of the Mira 

The analysis of the bolometric properties of the explosion re- 
quires a determination of the reddening independent of the XR 
analysis. Using the colors of the Mira variable, and an assigned 



Wavelength (A) 


EW (mA) 


E(B - V) 


5404 


13 


0.34 


5780 


234 


0.49 


6162 


27 


0.37 


6379 


29 


0.37 


6413 


54 


0.64 


6623 


109 


0.47 


6843 


12 


0.44 



Day 


kT (BB) 


kT (mekal) 


N„ 




after TO 


keV 


keV 


10^2 




<30 
30-60 
>60 


0.042+0.003 

064+""°2 


1 QQ+0.13 

0.80t!? 


0.62+0.09 




Day 


Obs flux 


Unabs flux 


Obs flux 


Unabs 


after TO 


total 


total 


hotter 


hotter 


<30 


6.17E-12 


2.44E-8 


5.81E-12 


1.09E-11 


30-60 


6.10E-12 


6.99e-ll 


5.59E-12 


1.17E-11 


>60 


1.34E-12 


6.14E-12 


1.34E-12 


6.14E-12 



Although there are no high resolution ultraviolet (UV) ob- 
servations of V407 Cyg with which to independently check the 
interstellar medium along the line of sight, the O star HD 202347 
lies at a close Galactic longitude and latitude. We show in Fig. 
3 a comparison of the V407 Cyg Na I D line with interstellar S 
II 1259 A resonance Hne in the STIS spectrum HD 202347, a Bl 
V star in NGC 7039. This is the closest OB star in Galactic lon- 
gitude and latitude in the currently available UV archives The 
velocity also agrees with that of the nearby (line of sight) B8 II 
star HD 199206, v^ad = -21.0 km s'' from Wilson (1953) and 
DDO Radial Velocity Catalog. This clearly identifies the low ve- 
locity component in the V407 Cyg spectrum as interstellar and 
also shows that the higher velocity components are from the star, 
especially the wind and circumstellar material. 

Several absorption lines, visible in the first NOT spectra, 
have been used to estimate the continuum contribution of the 
shock. These are from low-ionization resonance lines, known to 
be present in the quiescent spectrum and detected in the spectra 
of M6-M8 III stars. We show in Fig. 4 the comparison between 
the last NOT spectrum and stars from the ESO archive. The best 
match is to either M7 III or M8 III. The reported M6 III spectral 
classification of Munari et al. (1990) was obtained at a differ- 
ent phase (0.''1, around maximum light) of the Mira while the 
2010 outburst occurred at 0.''6 according to the Munari et al. 
ephemeris. 



3.3. Low-ionization species 

The O I emission features remained strong throughout the out- 
burst but the line profiles show a complex structure. In the first 
NOT observation, the [O I] 6300, 6364A lines displayed very 
narrow central emission peaks, FWHM =12 km s i, with low 



^ URL: http://arch ive.stsci.edu/| 
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Fig. 5. Variation of the [O I] 6300A line profile from the Ondrejov (11.4 km s ^px ^) and NOT (1.1 km s ^px ^) data sets during 
the observing period. The date is a ID-2450000; the peak of the XR emission, around day 30, corresponds to JD 2455296. 



05G301010, S II 1259.518, HD202347 vs V407 Cyg No I 5889-950 




V.loclly (km/s) 



Fig. 3. Comparison of the Na I Dl profile of V407 Cyg (dash) 
with the S 11 1257 A resonance line in HD 202347 (solid, archival 
HST/STIS G140H spectrum, program 05G3), along a close hne 
of sight. The high (negative) radial velocity components are from 
the Mira wind. The V407 Cyg spectrum, using the data from the 
first NOT observation, has been scaled for display; flux units for 
the HS r/STlS spectrum are erg s^^cm^^A"^ 
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Fig. 4. Comparison of the Ondrejov spectrum from MJD 55358 
with spectra for M6 III, M7 III, and M8 III stars from the ESO 
archive. 



level, extended wings. The core emission remained nearly invari- 
ant while the broad wings increased throughout the observing 
interval to a maximum velocity of 200 km s"^ This is shown 
in Fig. 5 from the NOT and Ondrejov sequences. The FWZI 
is nearly the same as the broad emission on the Na 1 D dou- 
blet, beginning at around the time of the break in the XR and 



UV emission seen in Fig. 1. The O I 8446A Une is fluoresced 
by LyP through the O I 1027 A resonance line (Bowen pump- 
ing), while [O I] 6300A is produced by recombination and cas- 
cade from absorption of the O I 1302A ground state transition, 
see Shore & Wahlgren (2010) for further discussion. The 8446A 
triplet narrowed during the interval from MJD 55279 (vrad,max=- 
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Fig. 6. Variation of the [O I] 6300A (open squares) and O I 
8446A (open diamonds) equivalent widths during the outburst 
based on the Ondfejov spectra. 



900 km s-i) to MJD 55315 (-400 km s"') after which it was un- 
observable in the Ondfejov spectra, being overwhelmed by pho- 
tospheric absorption bands. We display in Fig. 6 the variations 
of the equivalent width for [O 1] 6300A and O I 8446A. Note 
that the peak of the O I 8446A equivalent width is coincident 
with the XR peak and the increase of the [O I] 6300 coincides 
with the break in rate of decline of the UVOT fluxes and the 
6700A continuum. The maximum velocity observed in the O I 
8446A Une is about 500 km s ' . This suggests that the H I LyyS 
pumped transition, 1 1025 A, was as broad and was at its max- 
imum at the XR peak or that the combined effects of ionization 
of the environmental Ly/3 and the decrease of the intensity of the 
1 1025 A combined to produce the observed variation. No lines 
were detected for O II in the NOT spectra. 

The [O III] 4959, 5007A forbidden lines showed identical 
profiles and were initially similar to | O 111] 4363 A. The 4363 A 
Une, however, never developed the double peak structure shown 
in Fig. 7, displaying only the broad component. This indicates 
two separate contributors to the nebular lines, one from material 
at comparatively large distance that has not been swept up by 
the shock with a mean density of lO^cm"^ for 10 - 20kK, 
and the other, broader component with a density about a factor 
of 10 higher in the velocity interval consistent with impact of the 
ejecta. 

The [N II] lines, in contrast, show the two narrow compo- 
nents on aU three broad profiles (i.e. 5754, 6548, 6583A) al- 
though the latter two are invisible during the first three NOT ob- 
servations due to overwhelming blending with Ha (see Fig. 7). 
They were, however, detectable in the Ondfejov spectrum taken 
on MJD 55308, at the peak of the XR emission. Their intensi- 
ties increased relative to Ha by about a factor of 5 between MJD 
55315 and MJD 55334 and the [N II] (6548+6583)/Ha' ratio re- 
mained nearly constant after MJD 55352. The equivalent width 
ratio of [Ni II] 6667 A to He 1 6678 A, which is independent of the 
extinction and absolute flux calibration, shows a behavior simi- 
lar to that of the UVOT and 6700A continuum. The ratio steadily 
decreases throughout the outburst. 

The [Ar III] 7135.79A Une displayed the same profile as 
[N II] and [O HI] 4959, 5007 A, with two separate contributing 
regions. The photoionized region produced the same compara- 
tively narrow emission peaks and the [O I] 6300A Une has its 



Fig. 7. Comparative development of low ionization species Une 
profiles: [N II] 5754A (left column), [O I] 6300A (center), and 
[O III] 5007A (right) from NOT spectra on the three indicated 
dates during the outburst. 



miiximum at the local minimum between the narrow line peaks 
at -67 and -29 km s~' . The [Ar III] line, however, also developed 
a broad wing similar to [Ca V] (see below) extending from -200 
to -1-300 km s~' that contributed an approximately equal flux as 
the narrow Unes. In contrast, [Ar IV] 471 lA showed a similar 
profile to that of [Ca V] and He II but its peak is shifted in all 
spectra to H-110 km s"' (see below, sections 3.4 and 3.6). 

The Ca II IR triplet (8498, 8542, 8662 A ) lines all showed 
broad wings from the earliest Ondfejov spectra, extending to 200 
km s"' (Fig. 8). There were several narrower emission peaks in 
all profiles at the same velocity as the Na 1 D absorption com- 
ponents. In contrast, the [Ca II] 7291, 7323 A lines remained ex- 
tremely narrow, compatible with the [O I] 63 00 A core (12 km 
s"'), and never developed the broad wings that are seen on the Ca 
II 3968A resonance line or the fluoresced O I lines (Fig. 9). The 
Ca II Unes showed two very difl'erent developments during the 
outburst. The 8498, 8542, 8662A fi-iplet displayed a broad profile 
in the first observations that systematically narrowed throughout 
the observing interval, similar to O I 8446A. In contrast, the Ca 
II 7291, 7323A doublet remained narrow, «10 km s'^ (FWHM), 
and centered on the Mira, and never displayed broad wings. 

The Na I D lines displayed a broad component on which 
were superimposed the narrow absorption features (both ambi- 
ent and interstellar, see below), extending to 150 km s"'. The 
broad component increased steadily in strength throughout the 
observing interval, with the greatest change occurring after XR 
maximum. The centroid of the broad line appears to be shifted 
by about -16 km s ' relative to the Mira but this may be an ar- 
tifact of the fitting of a symmetric function (Gaussian) to what 
may be an intrinsically asymmetric profile of the same form as 
He II or [Ca V]. In Fig. 10 we show the Gaussian profiles for the 
two indicated epochs. The gallery of line profiles in Fig. 24 (see 
sec. 3.7) shows the changes in the broad component width (along 
with the narrow circumstellar absorption and emission features). 

The Mg 1 4571 A profile consisted of two components, like O 
I and [O I] : a narrow feature coincident with the Mira rest frame 
and broad wings extending from -200 to H-lOO km s~' in the 
first NOT spectra. The Une never displayed a P Cyg component 
or any of the absorption features observed on the Na 1 D lines. 
It strengthened during the outburst while remaining constant in 



7 



S. N. Shore et al.: The 2010 outburst of the recurrent nova V407 Cyg. I 




-200 200 

vrod (km/s) Co II 7291 



400 600 



Fig. 8. Top: Variation of the main line of the infrared Ca II triplet 

8542A; solid: MJD 55279, dash: 55294, dot-dash: 55315, dot- 
dot-dash: 55362 (Ondfejov spectra). Bottom: Comparison with 
Ca n 7291; soUd: 55289, dash: 55314, dot-dash:55363 (NOT 
spectra). 




vrod (km/s) I 8446-4 



Fig. 9. Variation of Bowen fluorescence O I 8446 A Une for the 
same epochs as Fig. 8 (Ondfejov spectra). 
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Fig. 10. Gaussian line profiles fit to the Na I Dl profile of the 
NOT sequence on the dates indicated. The line FWHM is 110 
km s"^ for both and the centroid is displaced to -70 km s"^ 



velocity width. In Fig. 11 we show the variation of the [O I] 
6300A and Mg I fines corrected for interstellar reddening. 

The equivalent variations of the [Ni II] 6667A Une relative 
to the nearly adjacent He 1 singlet 6678A is shown Fig. 12. 
This comparison is unafi'ected by either reddening uncertainties 



60 80 
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Fig. 11. O 1 6300A equivalent width variations relative to Mg 1 
4571A with a reddening correction E{B - V)=0.5 appUed to the 
continuum. 
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Fig. 12. Equivalent width variation of He I 667 8 A singlet 
(squares) relative to [Ni II] 6667 A (triangles). 



or continuum variations. The He I line systematically decreased 
during the observing interval. In contrast, the [Ni 11] 6667A line 
increased in equivalent with showing a change in rate of increase 
at around MJD 55300 similar to that seen for [O I] 6300A (e.g. 
Fig. 6) and [S II] 6730A (see below. Fig. 26). 

We will present a detailed study of the variations of the Fe- 
peak neutral and singly-ionized species in the next paper (in 
preparation). However, for illustration and comparison with the 
other neutral and low ionization species, we show in Fig. 13 the 
P Cyg profiles observed on the narrow components of the Fe II 
RMT 42 lines observable in the NOT spectra. The absorption 
persisted throughout the outburst, even as the wings of the pro- 
files increased in maximum radial velocity to 200 km s~^ The 
profiles were asymmetric at first, with extended blueward wings, 
becoming more symmetric in later spectra. The broad wings are 
clearly from the shock while the narrow componwent - and the 
P Cyg absorption - is from the Mira wind and chromosphere. 

The Si II 6347A line displayed a broad profile throughout 
the observing interval with wings extending to nearly 400 km 
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Fig. 13. Low velocity P Cyg absorption on the principal lines of Fig. 14. Variation of Si 11 6347A line profile in the NOT se- 

Fe II RMT 42 4923A (soUd), 5018A (dot-dash), 5169A (dot) on quence. The first spectrum has a contamination from Ha that 

MJD 55286, before the peak of the XR emission. Note also the aflrects adjacent orders from the extreme saturation of the pro- 

blueward extension of the line profiles (not shown beyond -200 file. This changes the core from -150 to -1-150 km s~' but has no 



kms-^. 



efl'ect on the extended wings. 



s ' and no P Cyg structure at any time. The companion line at 
6371A was always blended with Fe 11 6369.37, 6371. 13A and 
was unmeasurable. The profile was very different from the Ca 11 
IR triplet, never showing either the multiple emission features 
or a strong asymmetry. The Si II 5041, 5056A lines showed 
broad profiles throughout the observing interval but these are 
composite, complex blends with [Fe II] 5039, 5043, 5060A. In 
the first NOT spectrum, however, they are dominated by Si II 
for which we obtain the central velocities are both -90+10 km 
s"' and FWHM of 300+18 km s ' from Gaussian fitting. The 
lines were nearly of equal strength, and vvad.max was ±400 km 
s"' . The subsequent Une profiles are too complex to permit such 
a characterization. Note that the displacement and line width is 
similar to that of the Na I D broad components. The line profiles 
are shown in Figs. 14 and 15. The comparison to He I 6678 A, 
the only strong singlet observed in the outburst, shows the strong 
similarity in the profile and evolution. 

The He I line profile variations are shown in Figs. 17 and 
18. Except for 6678A, no singlet line of neutral helium was ob- 
served. The profiles of the 6678A Une, and those of the triplet se- 
ries, evolved similarly. The first NOT spectrum displays a broad 
(nearly +500 km s"'), nearly symmetric profiles for the triplets 
and 6678A, but the subsequent development is diff'erent. The 
6678A profile narrows but remains symmetric while the triplet 
lines develop multiple narrow emission features and an asym- 
metric shape extending from -200 to +500 km s"' . The minimum 
at the peak of the lines corresponds precisely to the maximum of 
the narrow [O I] 6300A emission, as found for [N 11] and [O III]. 

The He I line profiles developed very diflferently than dur- 
ing the 2006 RS Oph outburst. In the series shown in Fig. 51 
of lijima (2009), in the first week of the RS Oph outburst the 
profile showed a nearly triangular shape with wings extending 
to >1500 km s"', accompanied by narrower emission peaks at 
-200 and +100 km s"', that had narrowed to 800 km s"' by 
day 35. The subsequent evolution is close to a power law with 
Vrad,max ~ -800(?/35)"' km s"' and similarly for the red wing. 
For V407 Cyg, the He I and He II line profiles varied similarly 
except that the 5875A line showed two narrow emission com- 
ponents at the Une peak centered on the velocity of the Mira. 







A i 










Fig. 15. Variation of Si II 6347A in the Ondfejov sequence. 
Compare this with Fig 17 for He 1 6678A. The profiles are unaf- 
fected by the contamination present in the NOT echelle spectra 
at any epoch. 
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Fig. 16. Variation of Si II 6347A equivalent width from the NOT 
spectra. The peak in the emission occurs just after XR maximum. 
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Fig. 17. Variation of He I 6678A in the Ondfejov sequence for Fig. 19. Maximum radial velocity variation for the Balmer lines 
comparison with the Si II 6347A profile evolution. from the NOT spectra. Square: Hor, triangle: H/3, diamond, Hy. 



The He II profile in RS Oph, from Fig. 50 in lijima (2009), dis- 
played multiple components throughout the outburst after day 
35 and its broad wings remained symmetric, in contrast to V407 
Cyg. As in V407 Cyg, the singlets were weak or absent in RS 
Oph 2006. Only narrow components are reported for any singlet, 
e.g. 4921A. The He I 6678, 7065A Unes in RS Oph decreased 
steadily in time while remaining nearly constant relative to Hfi, 
the same was observed for V407 Cyg. 



3.4. HI Balmer line profiles 

The gallery of NOT spectra for the Balmer lines is shown in 
Fig. 18. The maximum velocity reported for the first Balmer line 
profiles (Mar 14 = MJD 55269, C. Bull, private communication) 
was >2500 km s"^ for Hor and HjS. The first Ondfejov observa- 
tion, (Mar 25, MJD 55280) shows, instead, extended emission 
to 4000 km s"'. This measurement is consistent with the Mar 
14 spectrum since the highest velocity is only seen in very weak 
wings (the S/N ratio for the first Ondfejov Ha profile is about 60 
pix~^ at 4000 km s"' and > 1000 at line peak). The profiles show 
a low velocity P Cyg absorption that remains strong through- 
out the period covered by these observations for the Balmer a 
through 7 lines, at -54 km s~'. This is consistent with the rest 
velocity determined from the Li I 6707A and Ba II 4554A reso- 
nance lines. Weaker absorption lines of metallic species seen on 
the Balmer profiles, especially on Hj3, decreased systematically 
until MJD 55332, after which they are not detectable. The H^ P 
Cyg trough at low velocity decreased systematically and turned 
into emission around the date of the UVOT light-curve descent. 

Figure 19 shows the measured extreme radial velocities for 
the Ha, H/3 , and Hy profiles. Blending prevented measurement 
of the extreme of the red wing of Hy and the low signal-to-noise 
ratio at H5 made the measurements unreliable. The uncertainties 
are ±20 km s"^ based on the continuum fitting. These should be 
compared with the He 11 4686A measurements shown below in 
Fig. 26. 



3.5. Li 1 6707A and non-photospheric continuum emission 

The Li I (flQlA line was visible from the first spectrum for which 
we measured an equivalent width, Wx - 42+5mA. The evolution 
of the line indicates that it was not formed in the wind but, in- 
stead, is the photospheric absorption feature previously reported 
by Tatarnikova et al. (2003a,b), who reported Wa = 350 mA. 
Its radial velocity, -54 km s~', is consistent with the Mira being 
at minimum light. During the first two months, the Une steadily 
increased until, after MJD 55350, it stayed constant. For the last 
two NOT spectra (medium resolution; MJD 55370, 55393) we 
obtained an equivalent width of 350+10 mA. If we assume the 
absorption line strength was, in fact, constant throughout the 
outburst, then the change was due to a "contamination" by a 
masking continuum; thus the additional continuum accounted 
for x!90% of the light at 6700A during the earliest observation 
(MJD 55286) and is consistent with the observed decrease of 
3 mag in R during the interval. This is consistent with the si- 
multaneous increase in the atomic and molecular absorption line 
spectra longward of 6700A, with the gradual approach to the 
pre-outburst spectrum shown in Munari et al. (1990). The mea- 
surements are shown in Fig. 20. As we will discuss in sec. 4, 
we interpret this continuum as probably being due to thermal 
bremsstrahlung emission from the shock. 

3.6. High-ionization species 

Shortly after the detection of y-ray emission by Fermi, Munari et 
al. (2010) reported the presence of high-ionization species in the 
optical spectrum. Our NOT sequence, covering the same period 
at higher resolution, allows us to confirm several of the identifi- 
cations. The [Fe X] 6375A line was present on the first Ondfejov 
spectrum (Mar 25) and continued to increase in strength through 
May. After the XR peak, it rapidly declined while maintaining 
the same line profile, almost identical to He II 4686A. The line, 
initially nearly symmetric with maximum velocity of ~500 km 
s"', became progressively more asymmetric in the same man- 
ner as the heUum lines and virtually disappeared by Jun 2. The 
[Ca V] Unes, in contrast, remained visible throughout the ob- 
serving interval and increased in strength through the XR peak, 
after which they decUned less rapidly than [Fe X]. For [Ar X] 
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Fig. 18. Line profiles for H I Balmer Ha, n/3, Hy, and U6; He H 4686A; He I 6678A (singlet) and 5875A (triplet) for all NOT 
observations. The abscissa is uniformly radial velocity in km s~^ the individual lines and MJD are indicated in the heading of each 
panel. The ordinate is flux in arbitrary units. 
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Fig. 20. Li I 6707-derived continuum contribution during the 
outburst, based on NOT spectra. 



5535A, the detection is ambiguous because of the presence of a 
set of strong permitted Fe-peak lines superimposed on a broad 
but weak emission. The Une is Ukely present in the first NOT 
spectrum. 

Broad emission was first observed as early as Mar 25. 
Munari et al. (2010) report the first detection of [Fe X] on Mar 



17, only one week after outburst, with a FWHM of about 600 
km s"'. This contrasts with lijima (2009) who reports, for the 
2006 outburst of RS Oph, that [Fe XJ and [Fe XIV] appeared 
only after day 35. [Fe X] 6375 A showed a nearly symmetric 
profile before and at XR maximum, while [Ar IV] 4711 A dis- 
played a He Il-like profile with an extended blue wing with a 
blueward maximum Vrad=-150 km s"' in the stellar reference 
frame, [Fe X] 6375A extended from -400 to +600 km s ^ This 
difference in the profiles persisted, the maximum blueward ex- 
tension of [Fe X] was always greater than [Ar IV]. The [Ca V] 
5309, 6086A lines were strong and showed a He Il-like profile. 
These are compared in Figs. 21 and 24. [Ne IV] 4714, 4720A 
and [Fe XIV] 5303A were absent. The feature observed near 
[Ar X] 5533A is more likely an Fe II permitted line, based on its 
profile in the first NOT spectrum and no other [Ar X] line was 
detected. 

Figure 22 shows the [Fe Vll] line evolution in the NOT se- 
quence and in Fig. 23, the last three epochs are compared to a 
simple model. The profile was computed using a spherical ex- 
pansion Monte Carlo method described in Shore et al. (1994). 
The maximum velocity was fixed to 400 km s"', consistent with 
the He 11 line blue wing, and a ballistic velocity law was used. 
The only adjustment required is to shift the line centroid to -1-70 
km s"' . The best agreement was achieved with a fractional thick- 
ness, AR/R = 0.6 + 0.2. The model reproduces the profile best 
in the two earliest stages, during the initial deceleration of the 
ejecta, while the last shows an asymmetry on the red wing, likely 
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Fig. 21. The changes in line profiles for He 11 4686 (top), |Ca 
V] 5308A (middle), and [Ca V] 6087A (bottom) on the last 
two NOT observations. The similarity of the two sets of pro- 
files demonstrates that the asymmetry is a dynamical effect and 
not an optical depth efl'ect (see text). 




600 



Fig. 22. Time development of [Fe Vll] 5720A in the NOT spec- 
tra, from bottom to top in the same sequence shown in Fig. 24. 



due to a change in the emission measure of the lower density gas 
with respect to that of the approaching side of the shock in the 
denser wind of the Mira. 



3.6.1. O VI Raman line profiles 

The broad emission features 6825A and 7080A present in many 
symbiotics are Raman scattering lines of O VI, formed by the 
conversion to long wavelength of O VI 1032, 1037 A doublet 
photons by the nearby LyyS 1025A fine (Schmid et al. 1996, 
Schmid et al. 1999). Although the UV flux was several orders 
of magnitude greater than a normal symbiotic star and the ion- 
ization was sufficient to produce the O VI FUV doublet, no O VI 
Raman emission was detected at either spectral region during the 
time interval covered by these observations. This contrasts with 
the RS Oph 2006 outburst during which the 6825A line was de- 
tected (Brandi et al. 2009, lijima 2009) in the early stages after 
outburst (before day 60). Brandi et al. note that the line was ab- 
sent in quiescence. Our non-detection may be due to the orbital 
phase at which the eruption occurred. If the WD was at supe- 
rior conjunction, the part of the Mira wind in which the O VI 
conversion occurred may have been occulted. In addition, unlike 
RS Oph 2006, at the time when recombination yielded sufficient 




Fig. 23. The observed [Fe Vll] 5720A line profiles (from MJD 
55363, 55370, 5593) compared with the simulation described in 
the text. 



colunon densities in LyjS to produce the Raman scattering, after 
Jun 3, the long wavelength region of the optical spectrum was 
dominated by the photospheric spectrum of the Mira. 



3.7. Emission line fiaring: Na i 5889, 5895A and [S 11} 6716, 
6730A 

The absorption line evolution of the Na I D lines has been 
discussed recently for the circumsteUar envirormient of sev- 
eral supernovae (Patat et al. 2007, Simon et al. 2009). We ob- 
served a systematic evolution of the absorption of several com- 
ponents during the post-outburst interval. In the first spectrum 
(MJD 55287), there is a very weak feature at -95 km s"^ with 
Wa - 5 + 1mA. Other absorption components were at -63, - 
59, and -54 km s"'. These latter remained virtually unchanged 
throughout the entire sequence and we checked the variations 
relative to the interstellar absorption feature at -23 km s"'. No 
variations greater than 10% were measured. Instead, the -95 km 
s"' feature systematically grew in equivalent width, reaching a 
maximum in our last NOT observation (MJD55393) of 120+2 
mA. At the same time, we noted the increase in an emission line 
at +6 km s"^ that lagged the increase in the absorption by ap- 
proximately two weeks. The sequence of line profiles is shown 
in Fig. 24. 

We propose that this is the rear side of the same shell that 
is responsible for the absorption at -95 km s~\ and that the de- 
lay measures the light-travel time across the structure. Assuming 
the measured velocity (in the absence of a dynamical model for 
the origin of the shell), the delay corresponds to a distance of 
approximately 4 x 10"^ pc (»900 AU) and a tentative time of 
origin of about one century ago. This is relatively close to the 
epoch of the 1936 "novalike" event reported by Hoffmeister that 
has been responsible for the denomination "nova" for V407 Cyg 
in the subsequent literature (a misreading of the catalog entry). 
However, it is well known that symbiotic stars undergo out- 
bursts, presumably thermonuclear runaways, during which time 
the degenerate mass accreter develops a strong wind. If this hap- 
pened in the extended outburst in the 1930s, then we may be 
seeing that ejection through the Na I variations. The timescale 
for the increase in the absorption line then suggests an electron 
density of about 6x10^ cm~^ at that distance. Without know- 
ing the thickness of the shell, we can nonetheless estimate the 
mass-loss rate using the column density, finding for this compo- 
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Fig. 24. Evolution of the Na 1 D doublet during the outburst 
(proceeds from bottom to top) from the NOT spectra. The spec- 
tra have been normalized for display. Note the appearance of 
the absorption on both components of the doublet at 5887.9 and 
5894.6A and the narrow emission at 5890.0 and 5896.7A. See 
text for discussion. 

nent M » 4 x lO^^'Mgyr ^ Since the 1936 event lasted nearly 
a decade, this amounts to a comparatively low mass of about 
10"*Mo, of which at least some part must be dust. 

The [S II] 6716, 6730A lines displayed a different profile to- 
ward the end of the sequence than the first, narrow profiles with 
a centroid velocity of -60 km s"^, hence very low velocity with 
respect to the star, spanning -150 to -1-50 km s~^ They remained 
centered near the stellar velocity and never displayed the red- 
ward extension of the higher ionization lines. A curious feature 
of their late time evolution is an apparent flare (see Fig. 26) af- 
ter MJD 55350 during which time the profile did not change. 
Flares of other low ionization species have been noted in lijima 
(2009) for the RS Oph 2006 event. This is more consistent with 
a change in the ionization state of the medium than an impact 
of the shock. In fact, it appears that the [S II] 6716, 6730 dou- 
blet was not formed in the shock but rather in the Mira, near the 
base of the wind. It indicates, however, the interval over which 
the other high ionization species are "contaminated" by environ- 
mental emission. In particular, the peak of the He II line is partly 
contributed by the low velocity broad emission and this accounts 
e.g. for the difference between it and the [Ca V] and [Fe X] pro- 
files. 



4. Discussion: Toward a model of the outburst and 
the binary system 

The 1985 and 2006 eruptions of RS Oph are the only symbiotic- 
like recurrent novae that have been followed throughout the 
event with modern, multiwavelength methods, and it is that sys- 
tem with which we will compare the 2010 V407 Cyg event. The 
optical Ught curves (see the supplementary material in Abdo et 
al. 2010) are almost identical although with different timescales, 
the V407 Cyg event was faster at all wavelengths. Although the 
spectroscopic development of the outburst was similar to that 
of the 1985 and 2006 events in RS Oph (Evans et al. 2008, 
lijima 2009) there were some important differences. The max- 
imum Balmer line velocity reached for V407 Cyg never ex- 
ceeded 3000 km s"' during the earliest stages of the expansion 
while those for RS Oph were well in excess of 4000 km s 




Fig. 25. Comparison between [S II] 6730A (dash), fCa V] 
6087A (dot), and He II 4686A (continuous) line profiles at the 
last NOT observation, MJD 55397. Note that the [S II] lacks the 
redshifted wing present on the higher excitation lines (see text 
for discussion). 
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Fig. 26. Variation of the [S II] 6730A equivalent with showing 
the shortUved "flare" event at about day 84 of the outburst. There 
was no profile change during this increase. Square: Ondfejov 
spectra, diamond: NOT spectra. Note the rapid increase in line 
strength approximately coincident with the peak and initial de- 
cline of the XR emission (see Fig. 1). 



The Balmer line profiles were also systematically different, be- 
ing more nearly broad Gaussians in V407 Cyg and nearly power 
laws in RS Oph (Anupama & Prabhu 1989, lijima 2009). While 
the line structure has been interpreted as individual emission re- 
gions, the V407 Cyg variations show that the low velocity P Cyg 
profile at around -90 km s"^ is absorption that changes to emis- 
sion late in the outburst (particularly evident on the Hi5 line). 
The persistence of the P Cyg absorption on both the Ha, H/3, 
Hy, and Fe II (permitted) lines at low velocity (around +10 km 
s~' in the system frame) is similar to that in the RS Oph 1985 
and 2006 outbursts. A difference is the switch to emission at the 
same wavelength for US at around the peak of the XR emission. 
The absorption was at the same velocity as the persistent narrow 
emission on the O I Unes, particularly the forbidden transitions. 

The extra continuum contribution derived using the Li I pro- 
file shows the same variation as the UVOT light curve. Scaling 
to the last date. Day 131 of the outburst, on which the contribu- 
tion is about 10% of the photospheric continuum, and using the 
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Fig. 27. Maximum positive (triangle) and negative (squares) ra- 
dial velocities for He II 4686A from NOT spectra. A correction 
of 50 km s"' has been applied to put the velocities in the rest 
frame of the red giant (see text for details). 



change in the Ba II 4554A line give a continuum consistent with 
bremsstrahlung assuming E(B - V)=0.5. On MJD 55363 (day 
104), for which we have a standard star calibration, the photo- 
spheric flux at 6705 A, near Li I, was 1 .6 x 10~'''erg s"^cm"-^A~^ . 
At maximum, in the first NOT spectrum, the extra continuum 
contributed a factor of 5=10 times higher (Fig. 20, bottom) at a 
time when the shock temperature from the XR mekal fit was 
» 2x lO^K, implying a velocity for the shock of ==800 km s ' for 
the interval from 30 to 60 days after TO. This is close to the mea- 
sured maximum velocity on the Balmer and He II 4686A lines 
at the epoch of the first two NOT observations. The shock could 
then have powered the high ionization emission with the optical 
continuum being contributed by cooler gas behind the expand- 
ing front, consistent with the emission from the lower ionization 
lines. This is also sufficient to produce the strong precursor that 
is required to account for the rapid appearance of the high ion- 
ization narrow components. The decrease in H I column density 
can be explained by the passage of the shock through the wind 
and the decrease in neutral absorber due to photoionization (see 
also Shore et al. 1996). 

The changes in the He I, He II, [Ca V], and Balmer line pro- 
files require an aspherical shock, the blueshifted side propagat- 
ing into the RG wind and the redshifted side displaying breakout 
of the ejecta. The terminal velocities can be fit with two power 
laws, the approaching side with an exponent of -0.84±0.05 and 
the receding side with an exponent of -0.31±.05. After about 
day 90, the blue side of the profile changes more slowly and 
is consistent with the continued expansion of a stalled shock. 
Assuming an adiabatic shock in self-similar expansion (Sedov- 
Taylor) with a power law, and an ambient density varying as 

with distance r from the Mira (constant mass loss, constant 
terminal velocity), the blueshifted side that is propagating into 
the Mira wind follows a velocity scaling as v ~ r^^'' while 



the redshifted side should vary as v 



-1/3 



. The highest ion- 



ization states (i.e. He^, Ca"^"*, Fe"^^, all of which have ionization 
potentials > 54 eV) display fine profiles that are nearly identi- 
cal among themselves but essentially different from those with 
lower ionization potentials (i.e. N^, O^, Ar^^, all < 54 eV). In 
particular, the highest ionization species, Fe^^, shows a third 



form for the line profile and also a systematic redshift. However, 
the blueward extension of the lines is asymptotically the same 
at around -200 km s~^ (i.e. «-150 km s"^ relative to the Mira). 
There are several possible explanations for the persistent blue- 
ward extension of the lines. One is the wake produced by the 
convergent shock on the rear side of the RG. This coUimated 
structure is a generic feature of the explosion and seen in both 
tiie Walder et al. (2008) and Orlando et al. (2009) numerical sim- 
ulations of the RS Oph outburst. The wake has a far longer life- 
time than the shock crossing time and is maintained by the hot 
post-shocked gas expanding from the cavity generated within 
the RG wind. The alternate, more prosaic, explanation is that 
the approaching side of the shock passed through the denser in- 
ner parts of the wind and continued to expand at constant ve- 
locity within the increasingly rarefied parts of the RG wind. The 
Balmer lines were more symmetric and consistent with emission 
from the entire shock front, while the highest ionization lines 
show the strongest redward extension and can be explained as 
the receding portion of the bubble (partly obscured by the Mira) 
from the low density outer wind. The lowest ionization species, 
especially O I and Na I, also displayed asymmetric profiles but 
with a less pronounced velocity ratio and lower velocities and 
can be accounted for as the recombintion in the cooling shock. 
The virtual absence of H I singlets relative to the triplets, and the 
broad lines observed on Na I, O I, and Si II in the post-shocked 
gas may indicate charge exchange effects in a thin shock. This is 
supported by the relatively late appearance of the neutral broad 
fines, at a stage after the asymmetric He II 4686A line appeared. 
Similar efl'ects have been discussed in the context of the Cygnus 
Loop and other older supernova remnants (e.g. Leutengger et al. 
2010). We will examine this in more detail in the next paper. 

The similarity of the outburst of V407 Cyg to that of RS 
Oph is striking, notwithstanding the diff'erences in the properties 
of the systems. In this case, modeUng performed for RS Oph 
(e.g. Walder et al. 2008) can also serve for interpreting V407 
Cyg. The explosion should have been non-spherical, although 
possibly less so than for the much shorter orbital period systems 
RS Oph and T CrB. Wind confinement toward the orbital plane 
may not be so developed in V407 Cyg as in the other two and 
it is plausible that the ejection was less collimated. Our simple 
line profile modeling of [Fe VII] supports this hypothesis. There 
is evidently neutral and weakly ionized material that has neither 
been highly ionized nor swept up in the shock. However, the dis- 
appearance of the narrow absorption lines on H^S and the changes 
in the high ionization species' profiles suggests that substantial 
acceleration has taken place in the ambient matter. Numerical 
simulations (Walder et al. 2008) displayed a near doubfing of 
the ejecta mass through mixing with the wind as the shock tra- 
verses this ambient medium within the first week. In the V407 
Cyg system this may be longer. 

As to the binary system itself, some quantitative constraints 
can be placed based on the similarity to RS Oph and the proper- 
ties of the RG. The Mira has a long pulsation period and based 
on the Munari et al. (1990) radius, 400 R©, it must not be fill- 
ing its Roche surface. If the strong Li I line is evidence of hot 
bottom burning, as argued by Tatarnikova et al. (2003b) (see 
also Herwig (2005) for a general review), then its mass must 
be greater than 4 Mq. Assuming that the WD is as massive as in 
RS Oph, around 1 .2 M©, then using the scaling from Eggleton 
(1983) for the ratio of the Roche radius, Rrl, to the semimajor 
axis, a, gives /?rl/« ~ 0.3 for a mass ratio of > 3. If Rrl > Ruka, 
then a > 6AU, consistent with the Fermi y-ray detection at 3 
days after optical outburst (Abdo et al. 2010) for a shock ve- 
locity of around 2500 km s ^ The period would therefore be 
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long, as suggested by Munari et al., (1990), Poib ^ 32 yrs. This 
is a far longer orbital period than either RS Oph (Brandi et al. 
2009) or T CrB, neither of which hosts a pulsating companion to 
the WD, and supports the conclusion that the degenerates in the 
symbiotic -like recurrent novae are the most massive among the 
symbiotic binaries Q 
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The referee pointed us to a study of RX Pup by Mikolajewska et al. 
(2002) that finds similarities between that system and both RS Oph and 
T CrB but with a longer period Mira variable (578 days) and a likely 
very long orbital period (>200 yrs). The line profile and photometric 
variations also show symbiotic nova-like outbursts similar to those ob- 
served in the 1930's in V407 Cyg. No explosive ejection has been de- 
tected from RX Pup but, in light of the noted similarities, it is worth 
monitoring. 



15 



